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ABSTRACT 

We report the discovery of a unique z = 6.027 galaxy, multiply imaged by the 
cluster Abell 383 and detected in new Hubble Space Telescope ACS and WFC3 imaging, 
as well as in Warm Spitzer observations. This galaxy was selected as a pair of i- 
dropouts; its suspected high redshift was confirmed by the measurement of a strong 
Lyman-a line in both images using Keck/DEIMOS. Combining Hubble and Spitzer 
photometry after correcting for contamination by line emission (estimated to be a 
small effect), we identify a strong Balmer break of 1.5 magnitudes. Taking into account 
the magnification factor of 11.4 ± 1.9 (2.65±0.17 mag) for the brightest image, the 
unlensed AB magnitude for the source is 27.2±0.05 in the H band, corresponding to 
a 0.4 L* galaxy, and 25.7±0.08 at 3.6 /xm. The UV slope is consistent with (3 ~ 2.0, 
and from the rest-frame UV continuum we measure a current star formation rate 
of 2.4±1.1 Mq/yi. The unlensed half-light radius is measured to be 300 pc, from 
which we deduce a star-forming surface density of ~10 M yr _1 kpc -2 . The Lyman-a 
emission is found to be extended over ~3"along the slit, corresponding to ~5 kpc in 
the source plane. This can be explained by the presence of a much larger envelope of 
neutral hydrogen around the star-forming region. Finally, fitting the spectral energy 
distribution using 7 photometric data points with simple SED models, we derive the 
following properties: very little reddening, an inferred stellar mass of M*=6 10 9 M , 
and an inferred age of ~800 Myrs (corresponding to a redshift of formation of ~ 18). 
The star-formation rate of this object was likely much stronger in the past than at 
the time of observation, suggesting that we may be missing a fraction of galaxies at 
z*~6 which have already faded in rest-frame UV wavelengths. 
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1 INTRODUCTION 

One of the most important topics of modern cosmology 
is the nature of the sources that reionized the intergalac- 
tic medium in the early Universe. A lot of progress has 
been made in characterizing the bright end of the popu- 
lation of high-redshift galaxies through the selection of opti- 
cal dropouts using both space- and grou nd-based data (e.g. 
iBouwens et al.ll2006l : TMcLure et al.ll2010l ). In particular, the 
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luminosity function of galaxies at z ~ 6 is well constrained 
down to 0.1 L* thanks to very deep observations of fields 
such as GOODS and the Hubble Ultr a Deep Field (UDF) , 
which provided over 500 i-dropouts (|Bouwens et al.1 120061 . 
120071 ). 

Solely based on photometric selection, the observed 
trend implies highly significant evolution at the bright end 
of the luminosity function (> L*) between z ~ 4 and z ~ 7, 
with a large decrease in the number of luminous sources. If 
this trend extends to higher redshifts, there might not be a 
sufficient number of ionizing photons to explain the observed 
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change of the intergalactic medium. One possible solution 
to this problem would be the presence of a high number of 
low-mass/low-luminosity sources at these redshifts, which 
would contribute to the bulk of the ionizing radiation. How- 
ever, only a small fraction of the currently known dropouts 
at z ~ 6 have luminosities of 0.1-0.5L* (L* corresponds to 
~26.3 AB at z ~ 6), and even fe wer are spectroscopically 
confirmed (|Stark et al.ll2010l . l201ll ). 

The newly installed Wide Field Camera 3 (WFC3) on- 
board the Hubble Space Telescope (HST) enables deep near- 
infrared imaging that is very fast compared to any other in- 
strument. So far most of the high-redshift work using WFC3 
has been performed on deep fields, the UDF in particu- 
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have been identified. However, the typical magnitudes of 
these objects (~27-28, and ~ 29 for the z ~ 10 candi- 
date) make ground-based spectroscopy challenging. Combin- 
ing WFC3 and the Spitzer Space Telescope (SST) at longer 
wavelengths allows us to study the star-formation history 
of these high-redshift galaxies, and to derive constraints on 
the age of the stellar population and the total stellar mass, 
thank s to the high sensit ivity of the first two IRAC channels 
(e.g., lEgami et ai]|2005l ). In particular, IStark et all (|2009h 
measured the masses of optical dropouts from z ~ 4 and 
z ~ 6, and found no significant evolution in the ratios of 
stellar mass to UV light. However, ages and stellar masses 
are best constrained in the brightest IRAC detections and 
quite uncertain for faint objects. 

A complementary and very powerful approach is to 
observe the high-redshift Universe through massive cluster 
lenses. Gravitational lensing offers the opportunity to per- 
form high signal-to-noise observations of galaxies that would 
otherwise be too faint and distant for quantitative study. Al- 
though the intrinsic lensed volume probed by any individual 
cluster observation is small compared to its unlensed equiva- 
lent, gravitational amplification by clusters likely represents 
a unique way to perform ground-based spectros copy of ob- 
jects that are intrinsically as faint as mag~ 30 ijEllis et al.l 
l200ll ; iKneib et alj|2004 IRichard et alll2008h . 

Combining the resolution and sensitivity of the Wide- 
Field Camera 3 (WFC3) and the Advanced Camera for Sur- 
veys (ACS) with the magnification afforded by clusters used 
as gravitational telescopes the recently approved Multi-Cycle 
Treasury "CLASH" cluster program on HST (PI: Postman; 
25 clusters; 524 orbits) offers a fresh opportunity to find 
and study high-redshift galaxies. Two thirds of the CLASH 
clusters have also been observed with SST/IRAC as part 
of our SST Warm Mission Exploratory Science program, 
"The IRAC Lensing Survey" (PI: Egami; 53 clusters; 593 
hours), which obtains deep (5 hours/band) images in the 
two shorter- wavelength IRAC channels (3.6 and 4.5/im). To- 
gether, these two large HST and SST observing programs fa- 
cilitate a detailed study of the physical properties of distant 
lensed galaxies. 

In this letter, we report the discovery of a high-redshift 
i-dropout in Abell 383 (ID 12065), which we confirmed to 
be at z = 6.027 based on Keck spectroscopy. We summarize 
the observations performed in Section 2 and present results 
of our photometric and spectroscopic analysis in Section 3. 
Conclusions are given in Section 4. Magnitudes are quoted 



in the AB system. At z — 6, 1" on the sky corresponds to 
a projected distance of 5.710 kpc in the adopted cosmology 
(On = 0.3, Ov = 0.7, h = 0.7). 



2 IMAGING AND SPECTROSCOPY 

HST images of the core of the massive galaxy cluster 
Abell 383 (a=02:48:03.332, $=-03:31:44.98,2 =0.187) have 
been obtained with ACS and WFC3/IR between November 
2010 and February 2011 in 10 filters from optical to near- 
infrared (F435W, F475W, F606W, F625W, F775W, F814W, 
F850LP, F110W, F125W, F160W) with exposure times of 
2000-3000 s per band. Individual ACS frames were corrected 
for charge transfer effici ency following the procedure de- 
scribed by Massevl (|2010h . adjusted to correct the most re- 
cent ACS images (more details will be given in Paraficz et 
al. in prep.). We measured the relative shifts between differ- 
ent epochs of observations in a give n filter and combined the 
ACS and WFC3 with multidrizzle (|Koekemoer et al.ll2002T l 
to obtain the final images. 

We then performed a search for i dropouts meeting the 
colour criteria F775W-F850LP> 1.3 l|Bouwens et all 12006ft 
in the strong-lensing region of the cluster. The envelope of 
the Brightest Cluster Galaxy (BCG) has been subtracted 
beforehand, using an elliptical surface-brightness model con- 
structed with the IRAF task el l ipse following the tech- 
nique presented in IRichard et al.l (I2008T I. This model pro- 
vides a good subtraction overall, except for residuals at the 
very centre that are due to the presence of a co mpact source 
and a dust lane, features already pointed out by Smith et all 
(|2005h . We used SExtractor l|Bertin fc Arnoutdll996h in 
double-image mode to obtain photometry across all HST 
bands. The F110W (most sensitive observation) was used 
as a detection image and the other bands as measurement 
images. Colours were measured between two bands inside 
the same 0.3" aperture and scaled to match the total mag- 
nitude in the F110W band (MAG_AUTO with an isophote 
of ~ 0.6" diameter). We used a 3-cr upper limit for the pho- 
tometry in case of non-detections. 

We identified two i-dropouts located 10" north (im- 
age 5.1: a=02:48:03.264, $=-03:31:34.77) and 25" south- 
east (image 5.2: a=02:48:04.600,$=-03:31:58.47) of the BCG 
(Fig. [TJ. They show a strong continuum break (F775W- 
F850LP > 1.9 mag.) and their colours agree within each 
other to better than the la error. The background gradient 
around image 5.1 is efficiently removed and the increased 
sky noise is estimated locally by SExtractor. The mea- 
surements of both images are listed in Table [1] Our well 
constrained lensing model for the mass distribution in this 
cluster (see Section[3J) predicts that both images 5.1 and 5.2 
are associated with the same source (referred to as source 
5), when assuming a redshift z > 5.5. 

On January 3, 2011 we used the DEep Imaging Multi- 
Object Spectrograph (DEIMOS) on the Keck II telescope 
to obtain an optical spectrum with a 1.2"-wide slit aligned 
to cover both images (Fig. [T]). Aiming to detect Lyman-Q 
emission from this source, we used the red-efficient 600ZD 
grating for three exposures of 1800 s. The reduction of the 
resulting data was performed with standard IRAF routines 
for bias subtraction, flat-fielding, sky subtraction and wave- 
length calibration. The standard star Feige 110 observed 
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Table 1. Summary of the measured photometry for both images. 
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Figure 1. Left: F814W-F110W-F160W colour image of the very core of Abell 383, together with known multiply-imaged systems 
(marked 1 to 5) and the critical line at z = 6. The identified i-dropout is seen as two images marked by yellow circles; the DEIMOS 
long-slit used in our spectroscopic follow-up is shown in white. Right: postage stamps of both images. HST and IRAC images are 3"and 
6"a side, respectively. 



during the same night was used for flux calibration. Seeing 
was l.l"and conditions photometric. 

A strong emission line is detected at 8545 A at the lo- 
cations of both images in the sky-subtracted spectrum (Fig. 
|2}. We measure an integrated flux of 2.2 ± 0.3 10~ 17 erg 
s _1 cm -2 in the extracted ID spectrum of the brightest im- 
age (5.1) and a spatial extent of ~3", much larger than the 
seeing of the observations. The line profile shows a possible 
asymmetry, but due to the close proximity of a skyline the 
signal-to-noise is too poor for an accurate measurement of 
the profile on the red side. However, all available evidence 
(including the lensing configuration and photometry) sup- 
port the identification of this line as Lyman-a at a redshift 
Zspoc = 6.027 ± 0.002 (measured at the peak of the line). 

In addition to the HST images, we used the SST/IRAC 
data of Abell 383 obtained at three different epochs (August 
2005, September 2009, March 2010). The first data set was 
obtained by programme 83 (PI: G. Rieke) for a total inte- 
gration time of 40 minutes/band. The second and third data 
sets were obtaind by programme 60034 (PI: E. Egami; "The 
IRAC Lensing Survey") and yielded a total 2-epoch inte- 
gration time of 5 hours/band. The Basic Calibrated Data 
(BCD) frames were mosaicked using MOPEX (MOsaicker 
and Point source Extractor) , a software package distributed 
by the Spitzer Science Center, using a pixel scale of 0'.'6 (half 
of the original instrumental pixel scale). 

Images 5.1 and 5.2 are detected in the two first IRAC 
channels (3.6^im and 4.5pm). We used a model of the BCG 
convolved with the IRAC point-spread function (PSF) to 



subtract its contribution and used a larger aperture (2") to 
measure the colour between the WFC3/F110W band and 
the IRAC bands. Using a l"aperture (more similar to the 
isophote used to infer the total magnitude) gives consistent 
colours but higher error bars. The WFC3 and IRAC images 
were convolved with an em pirical IRAC PSF and a tiny- 
TIM l|Hook fc Stoehrl 120081 ) model of the WFC3 PSF, re- 
spectively. We checked this photometry matching between 
HST and SST using cluster members as well as background 
sources of known spectroscopic redshift. Again, the photom- 
etry for image 5.1 (scaled to the total magnitude in F110W) 
is listed in Table [1] IRAC colours were not measured for im- 
age 5.2 since we were unable to remove the contamination 
from nearby non-elliptical cluster members. As the observed 
IRAC colours are entirely consistent between the 3 epochs, 
the BCG contribution to the sky noise is within the mea- 
sured errors. We adopt an average of the values from the 
two best epochs (2 and 3). 



3 PHYSICAL PROPERTIES 

We now estimate the strong-lensing magnification factors 
of images 5.1 and 5.2 based on a parametric model for the 
mass distribution of the c luster core. This model is con- 
structed using LENSTOOlQ l| Julio et al-l l^OO?), and based on 
a similar model constrained by 4 multiply imaged systems 

1 http://lamwws.oamp.fr/lenstool/ 
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Figure 2. 2D long-slit spectrum as obtained with 
Keck/DEIMOS, covering both images 5.1 and 5.2 and showing 
a clear emission line at 8545 A for both 5.1 and 5.2. The 
inset presents the extracted spectrum of image 5.1 around this 
wavelength. 

(re ferred to as system s 1 to 4, Fig. [J) pres e nted previously 
by ISmith et al.1 i|2003l ) and iRichardet all l|20ld) . Pseudo- 
isothermal elliptical mass distributions are used for both the 
cluster-scale component and individual galaxies. We include 
in our op timisation the r ecent spectroscopic redshifts ob- 
tained by iNewman et al.l (|201ll ) for systems 3 and 4. The 
location of the new pair of images identified in this letter 
(5.1 / 5.2) is in excellent agreement (within 0.5") with the 
predictions from this model at z = 6. We therefore use this 
new information as an additional constraint (system 5) to 
improve the mass model. From the family of best models 
resulting from the optimisation, we derive a magnification 
factor /ii — 11.4 ± 1.9 and [i2 — 7.3 ± 1.2 for images 5.1 
and 5.2, respectively. The ratio A*i//i2 is in perfect agree- 
ment with the F110W photometry. A third image (5.3) is 
predicted but is demagnified and located at the very centre 
of the BCG, making its detection extremely challenging, if 
not impossible. 

Taking into account these magnification factors, we de- 
rive intrinsic properties for the original background source 
5. Its unlensed F110W (J band) magnitude (27.4 AB) corre- 
sponds to an ~0.4L* galaxy at z = 6, which makes it more 
representative of the galaxy population at this redshift than 
bright er objects, such as the i-dropout identified in Abell 
1703 (|Richard et al.ll2009T l which corresponds to a ~3.6£* 
galaxy. A particularly interesting feature of source 5 is the 
large magnitude break (1.50±0.10 mag) between the H band 
(F160W) and the two first channels of IRAC, passbands that 
straddle the location of the rest-frame Balmer break. In or- 
der to interpret the implications of this break for the prop- 
erties of the underlying stellar populations, we performed a 
Spectral Energy Distribution (SED) fit of image 5.1 (high- 
est S/N image) with a library of spectr al templates follow- 
ing the precepts discussed in detail by IStark et al.l {2009). 
The main parameters we aim to derive are the total stel- 
lar mass (M*), the current star-formation rate (SFR) and 
the age of the stellar population (T). To this end, we re- 
moved the contribution of the measured Lyman-a line in 
the F814W and F850LP photometry, and then fit Chariot 
& Bruzual(2007, S. Chariot, private communication) stellar- 
population synthesis models to the observed SEDs. We con- 
sidered exponentially decaying star formation histories of 



the form SFR(t) ~ exp(— t/r) with e-folding times of r = 10, 
70, 100, 300, and 500 Myr, in addition to models with con- 
tinuous star formation (CSF). T is restricted to the range 
of 10 to 1000 Myr (the age of the universe at z ~ 6). We 
use a Salpeter initial mass func tion (IMF) (ISalpeterl ll955T ) 
and the dust extinction law of ICalzetti et al.l (|2000T ) with 
0.0 < E(B - V) < 1.0. Finally, we allow the metallicity to 
vary between solar (Zq) and 0.2 Zq, but note that metal- 
licity is degenerate with reddening and has a second-order 
effect on the SED. Absor ption by t he in tergalactic medium 
is accounted for following iMeiksinl (|2006T ) . 

We next consider the possibility that the strong break 
between the H band and the IRAC channels is due to con- 
tamination by nebular emission lines. Strong Ha, H/3 and 
[OIII] lines would all fall into the IRAC bands and could 
boost the continuum, mimi cking a stronger Balmer break 
jSchaerer fc de Barrosl I2010T ). We correct for this effect as 
follows: the SFR estimated from SED fitting is consistent 
with the observed Lyman-a flux, when assuming an aperture 
correction of a factor of 2 and (conservatively) 1/3 case B 
for Lyman-a. We predict observed fluxes /h„ = 1.46 10 -17 , 
f H /3 = 5.1 10" 18 and / [0 m] = 6.7 10" 18 (cgs units, both 
[OIII] lines average from COSMOS galaxies, Zoubian et al. 
2011 in prep.). We therefore lower the fluxes by 0.06 and 0.04 
magnitudes at 3.6 fim and 4.5/xm, respectively, thereby elim- 
inating contamination from unresolved line emission. This 
correction is negligible and we observe no significant change 
in the best-fit model using this updated photometry. 

The observed SED and the best- fitting template are pre- 
sented in Fig. [3] The best fit is achieved for a r = 500 Myr 
star-formation history. The full set of templates gives an 
age range of T=640-940 Myrs (within a A\ 2 = 3 from the 
best model) , depending on the exact star-formation history, 
corresponding to a redshift of formation of zj = 18 ± 4. 
More punctuated star-formation histories (r — 10 — 100 
Myrs) could produce similar Balmer breaks with ages of 
only T = 130 — 300 Myrs, but they produce much poorer 
fits to the data (A% 2 ~ 50, Fig. [3]). However, we acknowl- 
edge these are simple model s and age-dependent dust extinc- 
tion IjPoggianti et al.ll200lT ) as well as more complex star- 
formation histories (two-component models, with bursts on 
top of a continuous star-formation) could provide younger 
ages. This would affect the inferred ages more strongly 
than the inferred stellar masses, which are more indepen- 
dent of the assumed star-formation history. The most pre- 
cise measurement is the rest-frame V-band luminosity of 
2 - 2 -o'.42 lol ° L ©- The inferred values from the best fit model 
are a stellar mass of M*=6.3±i' 8 10 9 M , current SFR of 
3.2±1.1 M B /yr, and spe cific SFR (SSFR = SFR/M* , e.g. 
iBrinchmann et al.l 12004 ) of ~ 5 10~ 10 yr _1 indicate rela- 
tively vigorous sta r formation similar to the one found by 
lEgami et all ((2005). Note that the relative error in the mag- 
nification factor has been added in quadrature for all intrin- 
sic values. The rest-frame UV photometry, around 1500 A, 
is consistent with a model of constant AB magnitude, or 
equivalently with a slope of /3 ~ 2 (defined as fx oc A - ^). 
This is consistent with the very small amount of dust (E(B- 
V)< 0.06) in the best-fit model, and also with the result s 
found for z ~ 7 dropouts in the UDF (|Bouwens et al.ll2010h . 

Interestingly, image 5.1 of the galaxy is well re- 
solved in the WFC3/F110W image, as an ellipse of size 
~0.55" x0.3"based on well detected pixels. We measure a 
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Figure 3. Observed SED of image 5.1 (black data points) and 
model templates. The black curve represents the best fit to the 
photometry (r = 500 Myrs), whereas the blue/red curves show 
the best fit for younger but more punctuated star-formation his- 
tories (r = 10/100 Myrs respectively), giving a poorer fit to the 
photometry. 



half-light radius th = 0.21" with SExtractor, correspond- 
ing to a source-plane radius of 0.06", or 0.3 kpc after correct- 
ing for the PSF contribution. In comparison, the 3"extent of 
Lyman-a emission observed with DEIMOS, at least in the 
direction of the slit orientation, corresponds to a 5 kpc diam- 
eter in the source plane. This discrepancy consititutes strong 
evidence for the presence of a large envelope of neutral hy- 
drogen surroundin g the star-forming region responsible for 
the UV emission l|Steidel et al.|[201ll ). Of course, as these 
size measurements are given in UV rest-frame, it is quite 
likely that the source is imbedded in a lower surface- bright- 
ness halo, and the low reddening inferred would be coming 
from the UV-bright component. Nevertheless, a possible 3" 
diameter dusty component would have been resolved with 
IRAC. 



4 AN EVOLVED GALAXY IN A YOUNG 
UNIVERSE ? 

Taken at face value, the presence of a large optical break 
and an inferred mature stellar population (~800 Myrs) 
only 1 Gyr after the Big Bang is challenging. As pre- 
viously discussed, this value depends on the assumptions 
made in the standard m odels of star-formation history used. 
iMobasher et al.l (120051 ) reported a similarly strong Balmer 
break for a brighter and more massive galaxy at z ~ 6.5, 
which, however, still lacks spectroscopic confirmation de- 
spite numerous attempts. If we compare the intrinsic (i.e., 
demagnified) properties of this lensed galaxy with those of 
the 16 i-dropouts with spectroscopic redshifts analysed by 
lEvles et al l l|2007l ) and fit with similar SED models, this 
galaxy is 1 magnitude fainter in rest-frame UV (F850LP = 
27.4) but has a similar flux redward of the Balmer break 
(although with a much stronger detection thanks to its mag- 
nification by the cluster lens) and therefore the same typi- 
cal mass. The small size of the star-forming region implies 
a star-formation surface density of ~10 Mq yr" 1 kpc -2 , 
compar able to the most intense starburst activity observed 
locally l|Kennicuttl Il998l ) and allowing this object to drive 
strong winds. 

The detection of a large Balmer break, at such a high 



significance, is quite remarkable and demonstrates that sub- 
luminous systems (when selected in the UV) may not always 
be observed at a young stage of evolution. This is consistent 
with results o btained from IRAC stacking of faint z = 6 — 7 
field g alaxies ilstark et al.ll2009l ; lLabbe et ai1l2010| ), and 15% 
of the IStark et al.l "if2009) sample of bright i-dropouts show 
similar Balmer breaks. Thanks to the magnification pro- 
vided by gravitational lensing, we can measure this feature 
in faint individual objects. This result is important both for 
estimating the mass function at high redshift, and even more 
for understanding the process of reionisation, as it implies 
a significant contribution at z > 10 from the progenitors of 
these objects. 

Fro m our knowledge of t he UV luminosity function at 
2 ~ 6 l|Bouwens et al.ll2006 : ). we expect ~ 0.5 objects of 
this kind per CLASH cluster in the strong-lensing (multiply 
imaged) regime down to 0.2 L*. Deeper WFC3 and ACS ob- 
servations of strong-lensing clusters would efficiently make 
use of the Spitzer Warm mission program and allow us to 
measure the distribution of Balmer breaks in a larger sam- 
ple down to even fainter magnitudes, and thus to precisely 
determine the mass function. 
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